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Abstract 

The superconducting transition of a two-dimensional (2D) Josephson junction array exposed 
to weak magnetic helds has been studied experimentally. Resistance measurements reveal a 
superconducting-resistive phase boundary in serious disagreement with the theoretical and nu- 
merical expectations. Critical scaling analyses of the IV characteristics indicate contrary to the 
expectations that the superconducting-to-resistive transition in weak magnetic fields is associated 
with a melting transition of magnetic-field-induced vortices directly from a pinned-solid phase to a 
liquid phase. The expected depinning transition of vortices from a pinned-solid phase to an inter- 
mediate floating-solid phase was not observed. We discuss effects of the disorder-induced random 
pinning potential on phase transitions of vortices in a 2D Josephson junction array. 

PACS numbers: 74.25. Qt, 74.81.Fa, 74.25.Dw, 64.60.Cn 



When a two-dimensional (2D) Josephson junction array (JJA) is immersed in a perpen- 
dicular magnetic field, a finite density of vortices are induced in the array. The competition 
between the repulsive vortex- vortex interaction and the attractive periodic pinning potential 
due to the discreteness of the array results in interesting forms of superconducting-resistive 
phase boundaries as functions of a magnetic field". The superconducting-to-resistive tran- 
sition in strong magnetic fields or at high vortex densities has been discussed in terms of 
the melting of a vortex solid pinned to an underlying lattice . 2 ' 3 ' 4 ' 5 On the other hand, the 
superconducting-to-resistive transition in weak magnetic fields or at low vortex densities 
was found in simulation studies to be associated with the depinning transition of a pinned 
vortex solid, the appearance of which was predicted by Nelson and HalperinS in late 70s. 
Monte Carlo (MC) simulations of the lattice Coulomb gas model for vortices 7 showed that 
for a vortex density, the number of flux quanta per plaquette, / < 1/25 the pinned vortex 
solid transforms into a vortex liquid through two successive phase transitions: a depinning 
transition from a pinned solid to a floating solid at a lower temperature T p , approximately 
linear in /, and a melting transition from a floating solid to a liquid at a higher tem- 
perature T m ~ 0.045 J/ks, approximately independent of /, in which J is the Josephson 
coupling energy per junction. For / > 1/25, both transitions were found to appear at the 
same temperature increasing with /. These results have been confirmed semi quantitatively 
by different numerical or analytical studies of related models such as the XY model^, the 
resistively-shunted-junction model^, and the continuum model of a Coulomb gas in a peri- 
odic potential. While the existence of an intermediate hexatic phase just above the melting 
point has been also predicted by theory^i, no evidence was found in the numerical studies. 

Despite the agreement among the theoretical and numerical efforts, the experimental 
confirmation of the appearance of a depinning transition in the weak field limit has yet 
to be done. We, therefore, investigated experimentally the superconducting transition of a 
JJA in weak magnetic fields. The measurements of the resistance and the current-voltage 
(/V) characteristics of a square array represent a phase diagram in serious disagreement 
with the theoretical and numerical expectations. The expected depinning transition was 
not observed. We discuss effects of the disorder-induced random pinning potential on phase 
transitions of vortices in a 2D JJA. 

The experiments were performed on a square array of 400x600 Nb/Cu/Nb Josephson 
junctions. Cross-shaped 0.2-/im-thick Nb islands, fabricated by applying photolithography 
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and reactive-ion etching to the Nb/Cu double layer, were disposed on a 0.3-/im-thick Cu 
film periodically with a lattice constant of 13.7 /im, a junction width of 4 /im, and a junc- 
tion separation between adjacent sites of 1.5 /im. The photograph of Fig. C^a) shows the 
structure of the JJA sample used in this work. The variation of the junction separation in 
the sample was less than 0.1 /im. The standard four-probe technique was adopted for the 
measurements of the resistance and the IV characteristics. The current and voltage leads 
were attached to the Cu layer. The current leads were connected 3 mm away from the ends 
of the array to improve the uniformity of the injected current. Voltage leads were placed 100 
junctions away from the end edges of the array so that there exist 400x400 junctions be- 
tween the voltage leads. The sample voltage was measured by a transformer-coupled lock-in 
amplifier with a square-wave current at 23 Hz. The single-junction critical current i c at low 
temperatures can be obtained directly from the IV curve. At low temperatures, the slope 
dV/ dl of the curve has its maximum at a current equal to i c multiplied by the number of Nb 
islands perpendicular to the current direction. The i c and the junction coupling strength 
J(= hi c /2e) at high temperatures were determined by extrapolating the i c vs T data at 
low temperatures by the use of the de Gennes formula for a proximity-coupled junction 
in the dirty limit, i c (T) = z c (0)(l — T/T co ) 2 exp(— aT 1 / 2 ), where T co is the BCS transition 
temperature. We have also tried another formula developed by Zaikin and Zharkov^ and 
found both formulas generate practically the same z c 's for our sample at the temperatures of 
interest. Representing the superconducting transition temperatures of the sample in units 
of J/ks by using the z c 's determined with the de Gennes formula, they are 0.85 for / = 
and 0.42 for / = 1/2. These values are lower by 5-7% than the superconducting transition 
temperature (T c 's) found from many numerical studies, which are 0.89 for / = and 0.45 
for / = 1/2 13 . This demonstrates that i c {T) of the sample determined with the formulas 
is accurate within 10% of error. During the measurements, the temperature was controlled 
by a Lake Shore 340 temperature controller with fluctuations less than 1 mK. A solenoid 
generated external magnetic fields in the sample space where ambient magnetic fields were 
screened out by //-metal. The magnetic field or the vortex density was adjusted from the 
magnetoresistance measurements of the sample with 50-/xA excitation current. The presence 
of distinct resistance minima at fractional /'s in the magnetic field vs resistance curve made 
precise adjustment of the vortex density possible. 

Figure EJc) shows the sample resistance with 30-//A excitation current as a function of 
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/ at T = 3.59 K. The appearance of many pronounced higher-order dips indicates the 
good uniformity of the magnetic field over the sample. The magnetoresistance curve is 
impressively similar to the mean-field superconducting-resistive phase boundary numerically 
predicted in Ref. 1. The curve exposes, as numerically predicted, that the mean- field 
superconducting transition temperature T MF is a decreasing function of / in the low / 
region. As stated above, T c as an increasing function of / was observed in the same region 
in the MC simulations. Even though a mean-field phase diagram does not always resemble 
the real one, the /-dependence of T MF of our sample raises a question about validity of the 
phase diagram from the MC simulations. In order to examine the /-dependence of T c of 
a JJA in weak magnetic fields, we measured the sample resistance R as functions of the 
temperature T for six different vortex densities / = 1/50, 1/36, 1/25, 1/16, 1/10, and 1/8. 
Figure 121 shows some of the T vs R traces. The superconducting transition temperature T c 
determined from the T vs R curves is plotted against / in Fig. |3J The superconducting- 
resistive phase boundary differs from the theoretical and numerical expectations summarized 
above. T c decreases with / increased, contrary to the numerical observations. The T c at 
/ = 1/50 is ~0.32 in units of J/ks, which is ~10 times the numerical value of T p and ~7 
times that of T m . The T c 's at / = 1/50, 1/36, and 1/25 exceed even the melting temperatures 
of densely populated vortices for / = 2/5, 1/3, and l/5 4 i 5 i 14 . Such high T c 's appear hardly 
reconcilable with the proposed depinning transition, but possibly with the melting transition 
of a pinned vortex solid as for dense vortex systems. 

The gradual drop of the resistance of the sample upon approaching the transition indi- 
cates that the observed superconducting transitions are all continuous phase transitions. For 
a continuous superconducting transition, the critical scaling analysis of the IV characteris- 
tics may provide valuable information about the nature of the superconducting transition. 
Figure E| displays the IV characteristics of the sample for / = 1/50, 1/36, and 1/16. The 
log/ vs logV^ isotherms were obtained by averaging 15-240 measurements for each value of 
current. The derivative d(logV)/ti(logI), that is, the slope of the IV curves of Fig. 0] as 
functions of / are in Fig. The data at high currents are due to the single-junction effect 
and thus should be left out of consideration. The derivative plots show that for all the /'s 
studied, the IV curves have negative curvatures at low temperatures. At high temperatures, 
they have positive curvatures. The high-temperature and low-temperature curves are sep- 
arated by a straight line satisfying the power-law IV relation. Thermally activated vortex 
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motion would give rise to Ohmic IV characteristics in the low current limit and positive 
curvature in TV curves at intermediate currents. Negative curvature in the low-T curves 
implies vanishing resistance in the low current limit, indicating a superconducting state as 
the low-temperature state. The TV characteristics thus indicate that for the /'s, the JJA 
undergoes a continuous superconducting transition at the temperature where a straight IV 
curve appears. A detailed scaling theory^ suggests that for a continuous superconducting 
transition, the TV characteristics in 2D should scale as V/I\T — T c \ zu = £±(I /T\T — T c \ u ), 
where v and z are critical exponents and £± the scaling functions above and below T c . This 
scaling form becomes a simple power-law TV relation, V ~ J z+1 , at T = T c . Thus, one may 
find T c and the dynamic critical exponent z directly from the straight log/ vs logU isotherm. 
The correlation-length exponent v can be obtained from the scaling analysis of the IV data. 
Although the scaling analysis of IV data is an effective tool for studying the critical behav- 
iors of superconductors, special care should be taken not to extract any false information 
from the analysis. It has been recently shown for cuprate superconductors that the good 
data collapse can be achieved for a wide range of T c and critical exponents and sometimes 
does not even prove the existence of a phase transition.— It has been also shown that both 
the current noise^ and the finite-size effect may create Ohmic behavior at low currents 
even below the superconducting transition and lead to an underestimate of T c and incorrect 
v and z. The IV data in Fig. 0J however, do not suffer such complications. We avoided 
the current-noise problem by employing the phase-sensitive signal-detection method using 
a low-frequency square-wave current as described above. The effect of finite-size-induced 
free vortices was also found insignificant for JJA samples exposed to a magnetic field.— The 
IV curves in Fig. H] prove that the large flexibility in determining the critical exponents 
and temperature from TV data is not the case for the array sample. Unlike the TV data of 
cuprate superconductors, the IV curves of the array exhibit evident concavities below the 
transition. One can determine T c with uncertainty AT C « ±0.1 K from the curves. The T c 's 
determined from the straight IV curves agree within the uncertainties with those from the 
T vs R data. In the process of the scaling analysis, the distinct concavities of the IV curves 
further reduce the arbitrariness in determining T c to ±0.5-0.7 K. The TV data scaled on the 
basis of the scaling form with T c and z derived from the straight TV curves are shown in Fig. 
IH1 Each plot contains TV curves at 18-19 different temperatures. The scaling plots confirm 
that the JJA at the /'s studied undergoes a continuous superconducting phase transition 
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at the temperature where a straight IV curve appears. The insets of Fig. |H1 exhibit the 
values of T c , u, and z used to scale the data. The T c 's derived from the scaling analyses 
are in accordance within experimental errors with those from the resistance measurements. 
For all three /'s, the dynamic critical exponent z is less than 1 and the correlation- length 
exponent v is much larger than the 2D Ising value (uj = 1). The low-temperature IV curves 
can be fitted into the form V ~ Jexp[— (Jt//) m ] with \i = 0.7-1.1. The scaling behaviors 
of the IV characteristics are quite similar to those found for / = 2/5, 1/3, and 1/5&, for 
which the superconducting transition is understood in terms of melting of a pinned vortex 
solid driven by domain-wall excitations^*^. The large low-current voltage signals near the 
transition also appear to be compatible with the melting transition at T = T c . A floating 
solid is expected to be much less mobile near the transition and thus generate much lower 
voltage signal than melted vortices. Taking into consideration of the low vortex densities 
of our systems, the low-current voltage signals near the transition are as large as those of 
melted vortices for / = 2/5, 1/3, and 1/5 reported in Ref. 5. The large voltage signals above 
the transition, the scaling behaviors of the IV data, and the high T c 's, all indicate that the 
observed superconducting transitions in the weak magnetic fields are associated with the 
melting transition of a pinned vortex-solid to a vortex liquid without passing through an 
intermediate floating solid phase. 

The high melting transition temperatures and the absence of a depinning transition are 
in serious disagreement with the theoretical and numerical expectations. The disagreement 
in weak magnetic fields contrasts with the agreement in strong magnetic fields. For / = 2/5, 
the same sample exhibited a superconducting transition at about the same temperature as 
observed in simulations. It has been recently found in numerical studies of the 2D Coulomb 
gas^ that the floating vortex-solid phase may exist only in small arrays. This may explain 
the absence of a floating solid phase in our large array sample, but not the high melting 
temperatures. We therefore consider the disorder- induced random pinning potential as a 
possible cause for the disagreement in weak magnetic fields. The numerically obtained phase 
diagram 7 ' 10 intimates that the depinning and melting transitions may appear at the same 
temperature when the depinning temperature T p exceeds the melting temperature T m . A 
higher T p can be achieved with stronger pinning of vortices to the underlying lattice. In the 
process of fabricating a JJA sample, random variation of the junction coupling strength is 
inevitably introduced. The inevitable variation of junction coupling strength in a real array 
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may provide additional pinning and an enhanced T c for vortices dislocated from the periodic 
disposition. Such an effect of random bond disorder on T c should be much weaker in strong 
magnetic fields, consistent with the observations, because when densely populated vortices 
are dislocated from the periodic disposition, the increment of the repulsive vortex-vortex 
interaction energy effectively countervails the contribution of the random pinning potential. 
The variation of the coupling strength in our sample estimated from the / vs dV j dl curve in 
Fig. dh) is < ±15%. Yet it is not quite certain whether such an amount of variation of the 
coupling strength can raise T c so much. In order to ascertain whether extra random pinning 
can raise T c so effectively in weak magnetic fields, we measured resistively the T c of a system 
with artificially introduced random pinning. Figure 7 shows the T c as a function of / for a 
square JJA in which 14% of superconducting islands are randomly diluted. This site-diluted 
sample also contains the inevitable bond disorder, that is, the variation of junction coupling 
strength as much as the foregoing sample. We find in the figure that the random site disorder 
affects the T c of a JJA as presumed above. The T c in units of J/ks for / = 1/50-1/8 was 
elevated by the addition of site-disorder by 47-100% above that of the foregoing sample. For 
/ = 2/5, the change of T c was found less than 5%. The superconducting behavior of the 
site-diluted sample proves that the random pinning potential induced by disorder is quite 
effective in stabilizing the superconducting vortex-solid phase in weak magnetic fields. 

In summary, the measurements of the resistance and the IV characteristics disclosed a 
phase diagram of a JJA exposed to weak magnetic fields in sharp contrast to the theoretical 
and numerical expectations. Vortices in a square array at / = 1/50-1/8 were found to 
undergo a melting transition directly from a pinned-solid phase to a liquid phase at T c ~ 
0.15-0.32 J/Jzb, ~3-7 times the predicted melting temperature. The expected depinning 
transition was not observed. The experiments on the site-diluted array reveal that the 
disorder-induced random pinning can stabilize the superconducting pinned-solid phase in 
weak magnetic fields even at 10 times the melting temperature of a disorder-free system and 
is probably responsible for the experimental findings contrary to the theoretical expectations. 
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FIG. 1: (a) Photograph of the square Josephson junction array used in the experiment, (b) I vs 
dV/dl at T = 3.10 K. (c) Sample resistance with 30-fiA excitation current as a function of the 
vortex density / at T = 3.59 K. 

FIG. 2: Temperature dependences of the resistance at / = (a) 1/50, (b) 1/36, and (c) 1/16. The 
arrows denote where the superconducting transition occurs. The resistance was measured with 3 
\i A current injected into the sample. 

FIG. 3: Vortex-density dependence of T c determined from the resistance measurements. Note that 
T c is in units of J/ks- 

FIG. 4: The IV characteristics for (a) / = 1/50 at T = 3.600-4.500 K, (b) / = 1/36 at T = 
3.600-4.500 K, and (c) / = 1/16 at T = 3.500-4.400 K. The dashed lines are drawn to show the 
power law (V ~ I z+1 ) behavior at the superconducting temperature. 

FIG. 5: Current dependences of the slope of the IV curves in Fig. |1J Excessively fluctuating data 
at low currents are not shown. 

FIG. 6: Scaling plots of the IV curves. Each plot contains IV curves at 18-19 different tempera- 
tures. The insets show the values of T c , z, and v used to scale the data. 

FIG. 7: T c of the 14% site-diluted sample as a function of vortex density. 
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This figure "fig 1. JPG" is available in "JPG" format from: 



http://arXiv.org/ps/cond-mat/0506259vl 



This figure "fig2.JPG" is available in "JPG" format from: 



http://arXiv.org/ps/cond-mat/0506259vl 
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This figure "fig5.JPG" is available in "JPG" format from: 
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